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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 

in  Partial  Fulfillment  of  the  Requirements  for  the 

Degree  of  Doctor  of  Philosophy 

HEAT  TRANSFER  TO  WATER  DROPLETS  ON  A  FLAT  PLATE 
IN  THE  FILM  BOILING  REGIME 

by 

Kenneth  Joseph  Baumeister 

December,  1964 

Chairman:   Dr.  Robert  E.  Uhrig 

Major  Department:   Nuclear  Engineering 

The  mass  evaporation  rates  and  overall  heat  tranfer 
coefficients  are  determined  both  theoretically  and  experi- 
mentally for  water  droplets  which  are  supported  by  their 
own  superheated  vapor  over  a  flat  hot  plate. 

The  theoretical  and  experimental  mass  evaporation 
rates  are  found  to  agree  within  20  per  cent  over  a  drop- 
let volume  range  of  0.05  cc  to  1  cc  and  over  a  temperature 
range  of  600  F  to  1000  F.   In  this  parameter  range,  the 
mass  evaporation  rate  varies  from  0.001  to  0.01  (g/sec) , 
and  the  steam  gap  thickness  ranges  between  0.003  and 
0.008  in.   Both  the  gap  thickness  and  mass  evaporation 
rate  increase  for  increased  volume  and  temperature.   The 
overall  heat  transfer  coefficient  ranges  between  70 
(BTU/hr  ft2  F)  for  0.05  cc  droplets  and  40  for  1  cc  drop- 
lets in  the  temperature  reuige  considered.   Also,  the 

xii 


theoretical  analysis  yields  the  axial  and  radial  velocity 
distribution  under  the  droplet  and  a  velocity  correction 
factor  which  is  applied  to  Fourier's  equation  for  one- 
dimensional  steady  state  heat  conduction  across  the 
steam  gap. 

The  water  droplets  are  approximated  by  a  flat 
spheroidal  geometry  with  a  uniform  steam  gap  beneath  the 
droplet  and  a  saturated  steam  vapor  cover  on  the  top 
surface  of  the  droplet.   The  shape  of  the  droplet  and 
the  average  droplet  thickness  are  determined  analyti- 
cally.  The  analytical  results  compare  favorably  to 
experimental  measurements.   The  assumptions  are  made 
that  the  bottom  of  the  spheroid  is  at  the  saturation 
temperature  and  that  the  evaporation  takes  place  uni- 
formly beneath  the  spheroid.   The  flow  is  shown  to  have 
a  Reynolds  number  of  approximately  10;  consequently,  the 
flow  is  treated  as  incompressible  and  laminar  with 
negligible  energy  dissipation.   In  addition,  the  con- 
stant fluid  property  assumption  is  made,  and  because  of 
the  large  amount  of  time  required  for  the  evaporation 
of  the  droplet,  the  droplet  at  any  instant  is  assumed  to 
be  in  a  pseudo  steady  state  condition;  that  is,  the  flow 
is  approximated  by  a  steady  state  solution  at  any 

instant  of  time. 
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The  analytical  method  of  attack  is  to  solve  the 
momentum,  continuity,  cmd  energy  equations  simultane- 
ously.  The  partial  differential  momentum  and  continuity 
equations  are  reduced  to  ordinary  non-linear  differential 
equations  by  the  method  of  combination  of  variables. 
Possible  solutions  to  the  non-linear  equations  are  mapped 
by  means  of  an  analog  computer.   Then,  these  physically 
acceptable  solutions  are  combined  in  a  graphical  manner 
with  the  solutions  of  the  macroscopic  energy  equation, 
which  is  solved  explicitly,  to  yield  the  mass  evaporation 
rate  and  steam  gap  thickness  of  the  droplet  as  a  function 
of  droplet  size,  plate  temperature,  and  gravitational 
potential. 

The  effect  of  the  gravitational  potential  on  the 
mass  evaporation  rate  is  considered  in  detail  in  the 
theoretical  development.   A  reduction  in  the  gravita- 
tional potential  from  1  (earth)  to  0.16  (moon)  is  shown 
to  reduce  the  mass  evaporation  rate  by  approximately 
half. 
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CHAPTER  I 


INTRODUCTION 


The  object  of  the  study  presented  in  this 
dissertation  is  the  determination  of  the  overall  heat 
transfer  coefficient  from  a  heated  flat  plate  to  water 
droplets  which  are  supported  by  their  own  superheated 
vapor.   This  is  accomplished  by  a  theoretical  study  of 
the  momentum,  heat,  and  mass  transport  phenomena  associ- 
ated with  this  elemental  two  phase  flow  problem  amd  by 
an  experimental  verification  of  the  theory.   This 
analysis  differs  from  the  semi-empirical  and  dimensional 
approaches  used  in  the  past  in  that  the  analysis  is 
based  solely  on  the  solution  of  the  relevauit  governing 
equations  involved.   The  analysis  considers  a  wide  range 
of  surface  temperatures  and  volumes  of  water  droplets. 

Interest  in  this  subject  stems  from  the  rapid 
development  of  nuclear  reactors  as  used  in  power  and 
propulsion  systems.   In  particular,  there  has  been  a  new 
interest  awakened  in  the  general  problem  of  heat  trans- 
fer with  a  change  of  state,  such  as  occurs  in  the  two 
phase  flow  heat  transfer  in  a  boiling  water  reactor, 
nuclear  rocket,  or  in  the  more  fundamental  problem  under 
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consideration  in  this  dissertation.   Also,  an  Increased 
understanding  of  the  film  boiling  phenomenon,  such  as 
occurs  under  the  droplet,  is  of  importance  in  nuclear 
rocket  technology,  since  film  boiling  heat  transfer 
occurs  in  both  the  rocket  nozzle  and  core  reflectors. 
In  addition,  recent  works  by  Adadevoh,  Uyehara,  and 
Myers  (_1)   end  Borisheinsky,  Zamyatnin,  Kutateladze,  and 
Nemchinsky  (_2)  indicate  that  there  is  still  interest  in 
the  subject  of  droplet  vaporization  in  the  fields  of 
internal  combustion  engines  and  metallurgy. 

Investigations  on  this  subject  were  begun  as 
far  back  as  1756,  when  Leidenfrost  (3.)  first  described 
the  phenomena  of  film  boiling,  and  they  have  continued 
up  to  the  present  time  with  much  of  the  more  recent 
work  found  in  the  Russian  literature.   The  most  recent 
and  complete  works  are  by  Gottfried  (4^)  and  Borishansky 
C5). 

Gottfried  presents  both  a  dimensional  and 
semi-empirical  correlation  for  the  evaporation  of  small 
water  drops  on  a  flat  plate  in  the  film  boiling  regime. 


*The  underlined  numbers  in  parentheses  in  the  text  refer 
to  the  number  of  the  entry  in  the  List  of  References.   A 
statement  to  the  right  of  a  comma  within  the  parentheses 
will  give  the  location  within  the  source  to  which  the 
reference  is  made.   If  only  the  underlined  number  is 
given,  the  reference  is  to  the  source  in  general. 
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In  his  semi-empirical  approach,  the  water  droplet  is 
approximated  by  a  frustum  of  a  cone  in  which  the  upper 
and  lower  areas  are  varied  in  such  a  manner  so  as  to 
best  correlate  the  experimental  data.   From  the  upper 
area,  molecular  diffusion  is  assumed  to  take  place  in 
the  absence  of  thermal  convective  effects  or  a 
saturated  vapor  blanket.   From  the  lower  area,  uniform 
evaporation  into  a  superheated  vapor  is  assumed  with 
an  outward  flow  through  a  uniform  gap  beneath  the 
droplet. 

Borishansky  has  investigated  the  evaporation 
of  water  droplets  up  to  4  cc  in  volume  for  plate  tem- 
peratures of  527  F  and  662  F.   He  has  used  both 
dimensional  and  semi-empirical  techniques  in  correlat- 
ing his  data. 

On  the  bases  of  Borishansky ' s  experimental 
results  and  the  experiments  performed  in  conjunction 
with  this  dissertation,  the  general  problem  of  water 
droplet  evaporation  is  broken  down  into  the  following 
states,  which  are  governed  by  the  volume  of  the  drop- 
let:  small  spheroid,  flat  spheroid,  and  bubbly 
spheroid.   The  small  spheroid  state,  observed  at  0.05 
cc,  is  a  perturbation  of  a  purely  spherical  geometry 
by  the  action  of  gravity  working  against  the  forces  of 
capillary  tension  in  the  surface  of  the  droplet,  as 
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shown  in  Figure  1.   For  larger  size  droplets,  the 
perturbation  from  the  purely  spherical  state  increases. 
Finally,  as  the  liquid  volume  approaches  1  cc,  the 
droplet  reachers  the  flat  spheroidal  state  in  which  the 
thickness  of  the  spheroid  undergoes  little  change  with 
further  increase  in  its  volume.   For  volumes  greater 
than  approximately  1.5  cc,  the  vapor  formed  beneath  the 
droplet  tends  to  breeUc  intermittently  through  the  sur- 
face of  the  liquid  which  gives  rise  to  the  term  bubbly 
spheroid  state,  shown  pictorially  in  Figure  1. 

In  this  particular  study,  water  droplets  in  the 
volume  range  0.05  to  1  cc  are  analyzed.   In  this  volume 
remge,  an  analytical  model  based  on  a  flat  spheroidal 
geometry  reasonably  satisfies  the  physical  situation 
and  yet  still  has  simple  enough  boundary  conditions  to 
make  the  resulting  boundary  value  problem  tractable. 
The  shape  of  the  droplet  and  the  average  droplet  thick- 
ness, i,    are  determined  analytically.   The  analytical 
values  of  i   compare  favorably  to  the  measurements  taken 
by  Borishansky. 

The  analytical  method  of  attack  is  to  solve  the 
momentum,  continuity,  and  energy  equations  simultane- 
ously.  The  partial  differential  momentum  and  continuity 
equations  are  first  reduced  to  ordinary  non-linear  dif- 
ferential equations  by  the  method  of  combination  of 
variables.   A  difficulty  results  from  the  fact  that  the 
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6 
non-linear  equation  cannot  be  solved  in  closed  form. 
The  difficulty  is  further  compounded  because  the 
boundary  conditions  to  the  problem  are  unknowns.   The 
analytical  approach  to  the  problem  is  to  map  possible 
solutions  to  the  non-linear  flow  equations  by  means 
of  an  analog  computer  and  to  use  those  particular 
solutions  which  are  physically  acceptable.   Then,  the 
solutions  to  the  macroscopic  energy  equation,  which 
is  solved  explicitly,  are  combined  in  a  graphical 
manner  with  the  results  of  the  amalog  solutions  to 
yield  the  mass  evaporation  rate  of  the  droplet  as  a 
function  of  droplet  size  and  plate  temperature.   The 
overall  heat  transfer  coefficient  can  be  found 
directly  from  the  mass  evaporation  rate. 

The  theoretical  mass  evaporation  rates  are 
found  to  agree  within  20  per  cent  of  the  experimental 
values  over  a  temperature  range  of  600  F  to  1000  F, 
and  a  volume  range  of  0.05  to  1  cc.   The  analysis  also 
yields  the  axial  and  radial  velocity  distribution 
under  the  droplet,  the  steam  gap  thickness,  and  a 
velocity  correction  factor  which  is  applied  to  Fourier's 
equation  for  one-dimensional  steady  state  heat  conduc- 
tion across  the  gap. 

Also,  with  the  possibility  of  the  moon  being 
explored  in  the  next  decade  and  with  the  possibility 
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of  operating  with  two  phase  flow  in  a  low  gravitational 
field,  the  theoretical  analysis  takes  into  account  the 
effect  of  a  variation  of  the  gravitational  potential  on 
the  mass  evaporation  rate.   In  particular,  the  mass 
evaporation  rate  is  determined  for  the  physical  situa- 
tion where  the  surface  gravity  is  equivalent  to  that 
of  the  moon's  surface.   The  resulting  values  of  the 
mass  evaporation  rate  are  compared  to  the  values 
obtained  on  earth  for  a  similar  plate  temperature  and 
droplet  volume.   The  mass  evaporation  rate  is  approxi- 
mately half  that  found  on  the  earth  for  a  given  droplet 
volume  and  plate  temperature. 


CHAPTER  II 

METHOD  OF  ANALYSIS 

General  Approach 
Consider  the  flat  water  spheroid  shown  in 
Figure  2.   Heat  transfer  to  this  flat  spheroid  takes 
place  primarily  by  conduction  and  radiation  through 
the  superheated  film.   Heat  transfer  and  evaporation 
from  the  upper  surface  are  negligibly  small  in  com- 
parison to  that  beneath  the  droplet.   Kutateladze 
(6,  p.  376)  points  out  that  the  external  surface  of 
the  spheroid  is  covered  by  superheated  vapor  flowing 
from  beneath  the  spheroid.   This  vapor  cover  reduces 
the  energy  transport  from  the  upper  surface  to  a  near 
zero  value.   However,  even  in  the  assvuned  absence  of 
a  steam  cover,  both  the  free  convective  and  radiative 
heat  transfer,  and  free  convective  evaporation  are 
negligible  when  compared  to  that  occurring  beneath  the 
droplet.   The  free  convective  evaporation  was  esti- 
mated from  a  correlation  presented  by  Wado  (2) • 

Thus  the  problem  of  heat  transfer  to  a  flat 
water  spheroid  reduces  to  a  problem  commonly  termed 
mass  transfer  cooling  (8) .   Such  a  process  is 
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characterized  by  a  mass  flow  through  a  porous  surface 
(transpiration  cooling) ,  by  mass  released  from  a 
surface  through  evaporation  or  sublimation,  or  some 
chemical  reaction  (film  cooling,  ablation  cooling) . 
To  determine  the  velocity,  temperature,  and  mass 
distribution  for  the  above-mentioned  problems,  Hartnett 
and  Eckert  (9.)  point  out  that  in  general  it  is  neces- 
sary to  solve  the  continuity,  momentum,  energy,  auid 
diffusion  equations  simultcineously.   However,  for  the 
superheated  steam  region  under  the  water  droplet,  the 
mass  evaporating  from  the  surface  is  the  seune  as  that 
flowing  beneath  the  droplet;  thus,  the  diffusion  equa- 
tion is  identicallly  zero,  as  pointed  out  by  Grober, 
Erk,  and  Grigull  (10,  p.  416). 

However,  for  the  problem  of  evaporation  of  a 
flat  spheroid  resting  on  its  own  superheated  film,  an 
added  condition  of  static  equilibrium  is  required  for 
the  solution  to  the  problem.   The  pressure  forces  on 
the  bottom  of  the  droplet  must  be  sufficient  to  balance 
the  weight  of  the  droplet  minus  the  reactive  force. 
For  this  particular  problem,  the  reactive  forces  are 
negligible  compared  to  the  body  force  due  to  gravity, 
as  shown  in  Appendix  A.   Consequently,  the  reactive 
force  is  neglected  in  this  analysis. 
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Gottfried  (4)  and  Kutateladze  (6,  p.  377)  point 
out  that  the  flow  under  consideration  is  of  very  low 
velocity  and  is  well  within  the  laminar  range;  thus, 
the  flow  is  treated  as  incompressible  with  negligible 
energy  dissipation  due  to  friction.   In  addition, 
because  of  the  large  amount  of  time  required  for  the 
evaporation  of  the  spheroid,  the  droplet  at  any  instant 
is  assumed  to  be  in  a  pseudo  steady-state  condition; 
that  is,  the  flow  is  approximated  by  a  steady-state 
solution  at  any  instant  of  time.   Consequently,  for 
this  case  of  axisymmetrical  and  incompressible  laminar 
flow  with  negligible  dissipation  and  with  constant 
fluid  properties,  the  momentum,  continuity,  energy,  and 
static  equilibrium  equations  are  as  follows: 

Momentum: 


du  bn  9c    ^>P  /d^u    1  du    u  ^  6^u  ,   ,,, 

u  —  +w  —  = —   —  +  ^    + ^  I-  (1) 

Or  ^z     p  dr  \>)  2   r  dr        j-^   dz' 


u  —  +w  — —  +  v| +  _  __  +  J.       (2) 

dr  dz 


Continuity! 


^  +  ^  +  ^  -  0  .  (3) 

dr        r    dz 
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Energy X 


u  ^  +  „  |I  =  aV^'T  .  (4) 

dr  dz 


Static  Equilibrium: 


1° 


P(r,6)  iTfrdr   -  wf^l-   ,  (5) 

go 


The  boundary  conditions  for  the  above  equations  are: 


z  =  0     u-0     w-^0        T  =  Tp.     (6) 
2  =  5     u  =  0     w  =  w(6)     T  =  Tgat  •   <"7) 


The  assumptions  are  made  that  the  bottom  of  the 
spheroid  is  at  the  saturation  temperature,  and  that  the 
evaporation  takes  place  uniformly  beneath  the  spheroid. 
The  boundary  condition  on  the  axial  velocity  at  the 
upper  surface  is  an  unknown;  in  fact,  at  the  present 
time  the  gap  thickness  is  also  an  unknovm.   The  gap 
thickness,  O  ,  and  w((j)  ,  are  found  by  simultaneous  solu- 
tion of  the  above  equations. 

In  this  cuialysis,  the  determination  of  the 
evaporation  rate,  heat  transfer  coefficient,  aund  gap 
thickness  is  in  terms  of  the  volume  of  the  water  droplet, 
gravitational  potential,  and  temperature  of  the  heating 
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plate.   Consequently,  the  transport  properties  of 
viscosity  and  thermal  conductivity,  as  well  as  the 
specific  volume  of  the  steam,  are  expressed  in  terms 
of  the  steam  temperature.   For  the  range  of  tempera- 
ture under  consideration  in  this  dissertation,  the 
above  properties  are  represented  as  linear  functions 
of  temperature. 

Momentum  Equation 
The  logical  beginning  to  this  analysis  is  to 
solve  the  momentum  equation  since,  as  a  result  of  the 
constant  fluid  properties  assumption,  the  mutual 
interaction  between  the  equation  of  motion  amd  the 
energy  equation  ceases,  and  the  velocity  field  no 
longer  depends  on  temperature.   The  usual  approach  (9.) 
is  to  evaluate  the  properties  of  the  flow  field  at 
the  film  temperature,  as  defined  as: 

_  "^p  -^  "^sat  (8) 

f  "    2 

Therefore,  the  immediate  problem  is  to  solve  the 
continuity,  Equation  (3)  and  the  momentum.  Equations 
(1)  and  (2),  simultaneously.   These  equations  form  a 
set  of  partial  differential  equations  with  two  inde- 
pendent and  three  dependent  variables.   Use  of  the 
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method  of  "combination  of  variables"  at  this  time 
reduces  the  above-mentioned  equations  into  a  set  of 
ordinary  non-linear  differential  equations. 

Consider  for  a  moment  the  physical  situation. 
Defining  u  as  a  mean  radial  velocity, 


_  _  1    ^" 

u  -  -       u  dz  ,  (9) 


and  the  z  component  of  velocity  at  the  surface  of  the 
water  spheroid  as  w(6),  then  the  conservation  of  flow 
into  and  out  of  a  cylindrical  volume  of  radius,  r, 
under  the  flat  spheroid  results  in 


w(6)Tr  r^  =  u  2n  r6 


(10) 


However,  since  w(5)  is  an  assumed  constant  along  the 
bottom  surface  of  the  spheroid,  it  follows  that 


u  a  r 


(11) 


Also,  Gottfried  (4),  working  with  the  mean 
radial  velocity  as  defined  above  and  with  the  Navier- 
Stokes  equation  in  the  radial  direction  only,  shows 
that 
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The  above  relationships  indicate  that  a 
conbination  of  variables  of  the  form  used  in  the  prob- 
lem of  three  dimensional  axisymmetry  stagnation  flow 
^,     (11,  p.  83)  will  reduce  the  partial  differential 
^>^     equations  (1),  (2),  and  (3)  into  a  set  of  ordinary 

equations .   The  functional  forms  used  in  this  conver- 
sion are: 

w  =  -2  f(z),  (13) 

u  =  r  f  (z),  (14) 


P  =  1  a2  £. 


gc 


ro2  -  r2  +  F(2)   .  (15) 


Choosing  these  specific  functional  forms  for  w 
and  u  satisfies  identically  the  continuity  equation, 
as  substituting  Equations  (13)  cmd  (14)  directly  into 
the  continuity  equation  (3>  verifies. 

The  partial  derivatives  of  the  relationships 
(13),  (14),  and  (15)  used  in  the  combination  of  vari- 
ables substitution  are  now  listed  for  future  reference. 


(16) 
(17) 
(18) 


1^- 

f "  =  °  ■ 

dz    ^    ^ 

f  ' 

d2u 

A    2   -   ^   ^ 

^^  -   0 
dr 

dr2 
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^  "     ^  r-2  ^  -2  f"  '        (19) 

Oz 

^  =  -,2   "  r        ^^    1   "  2 

Substituting  the  above  relationships  into 
Equations  (1)  and  (2)   results  in  the  following  new 
forms  of  the  momentum  equations: 


f'2  -  2  ff"  =  a2  +  Vf 


(21) 


2  ff-  .  -ia2  F'  -  ^f..  .  (22) 


Thus,  the  original  partial  differential 
equations  are  by  means  of  the  above  substitutions  con- 
verted into  a  set  of  non-linear  ordinary  differential 
equations.   Although  these  equations  are  not  easy  to 
solve  in  the  strict  sense,  they  are  much  more  easily 
handled  than  the  original  partial  differential  equa- 
tions . 

Gottfried's  pressure  relationship   (Equation 
12)   is  generalized  to  include  z   variations  in  pressure, 
since  a  trivial  solution  to  the  momentum  equation  re- 
sults if  the  z  variation  in  pressure  is  neglected,  as 
shown  in  detail  in  Appendix  B.   The  boundary  conditions, 
Equations  (6)  and  (7),  now  take  on  the  following  forms: 
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z=0     f=0         f'=0    F=0    T=Tp,    (23) 
2  =  5    f  =  '  ^^^        f    =   0  T  =  Tgat  •  (24) 


The  boundary  condition  on  the  function,  F,  is 
chosen  arbitrarily  since  the  only  interest  is  in  the 
relative  variation  of  the  function  in  the  z  direction. 
The  constant,  a^,  accounts  for  the  absolute  magnitude 
of  tn  e  pressure. 

Equations  (21)  and  (22)  are  now  freed  of  the 
constants  a^  and  V  by  making  the  following  substitutions: 

C=C3^z     f(z)=C2  0(C)      P(z)  =  C3  !//(().   (25) 

Substituting  the  above  expressions  into  Equations  (21) 
and  (22) ,  the  coefficients  of  these  equations  become 
independent  of  a^  and  Z^  if 

c  =^/l  - .  (26) 


f(z)   =  Ya270(O  .  (27) 

F(z)   =  ^l//(0  .  (28) 
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Therefore,  it  follows  that 


d^  dz 


(29) 


(30) 


f  •   =  a-^d"    ,  (31) 

f"   =  -y-0'"  ,  (32) 


i; 


a  d'"  dz     a  yv^ 


(33) 


Substituting  the  above  relationships  into 
Equations  (21)  amd  (22) ,  the  momentum  equation  takes 
on  the  following  forms: 

(j)'^   -   2(p(p'  '  =  1  +  {)'  ".  (34) 

2(p(p'    =  -y/'  -  0"  .  (35) 

The  boundary  conditions  on  these  equations 
become 

(^   0  0-0        0'  =  0     yy=  0    T  =  T 

(36) 
(^  =  -  y?  /;    ^  =   -^(J)    r^.  =  n  T  =  T 

(37) 


^i^         2  Vai7 
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Complications  result  because  a  closed  solution 
to  the  above  non-linear  ordinary  differential  equation 
does  not  exist  and  because  6  euid  w(6)  are  unknowns. 
The  value  of  the  parameter,  6  ,  depends  on  the  simultan- 
eous solutions  of  all  the  governing  equations.   Conse- 
quently, the  method  of  solution  is  to  assume  many 
reasonable  values  of  *  and  to  solve  for  the  flow 
distribution  in  each  of  these  cases  that  satisfy  the 
static  equilibrium  condition.  Equation  (5) .   The  results 
are  plotted  and  compared  to  the  solutions  of  the  energy 
equation,  Equation  (4) ,  for  various  values  of  5  .   The 
intersection  of  the  results  of  the  momentum  cind  energy 
equations  represents  the  value  of  the  gap  thickness,  6.. 

The  solution  of  Equations  (34)  and  (35)  in  this 
particular  situation  is  most  easily  performed  by  means 
of  an  analog  computer.   However,  instead  of  assuming 
values  of  <5  ,  initial  values  of  4>  '  '  are  assumed.   The 
values  of  ^    and  w(5  J  are  determined  from  the  output  of 
the  computer. 

The  next  section  discusses  in  detail  the  analog 
solution  to  Equations  (34)  and  (35) . 
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Analog  Solution  of  Momentum  Equal- ion 
The  solution  of  ordinary  differential  equations 
by  use  of  a  differential  analog  computer  is  discussed 
widely  in  the  literature  (12 ) ,  (13 ) ,  and  (14) .   Briefly, 
the  analog  computer  is  capable  of  the  basic  mathematical 
operations  of  addition,  subtraction,  multiplication,  and 
integration.   The  variable  quantities  of  the  differen- 
tial equation  are  represented  by  voltages  which  may  be 
recorded  by  use  of  the  proper  recording  equipment. 

In  setting  up  the  momentum  equations  (34)  and 
(35)  for  solution  on  the  analog,  it  is  helpful  to 
rewrite  these  equations  in  the  form: 

0'"   =  0''  -2^p-'  -  1,  (38) 

^'      -      -   2(1)0'    -(])"  ,  (39) 

The  first  step  in  programing  the  above  equations 
for  the  analog  computer  is  to  change  the  variables  in 
the  above  equations  to  computer  variables  by  making  the 
following  transformations: 

C     =     Cf  T,  (40) 

^(0   =  <i>iT),  (41) 

V^^O   -  /?  1^^(70  -  (42) 
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The  derivatives  of  the  above  transformations 
necessary  for  substitution  into  the  auialog  equations 
(38)  and  (39)  are 

^.    =   £i£l   =   ±^   =  1:11  .=  1  •  (43) 


(44) 


d((,   dx 

1  d(^ 

K   d* 

K       . 

dT   dC 

0     (JT 

a   di 

0 

*•• 

= 

If 

1  I  t 

K 

3   $ 


iJ'       a   '*'  * 


(45) 


(46) 


Substituting  the  above  relationships  into 
Equations  (38)  and  (39)  results  in  the  analog  momentum 
equations 

V  =  <ai^  -  2<o*;   -^   ,  (47) 

i     =-2   S\    ;     -  JL-  ;  .  (48) 

The  boundary  conditions  on  the  above  equations 
take  on  the  form 


T=0  (l>  =  0  *  =  0    H'  =  0    T  =  Tp(49) 


^  (50) 
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The  symbols  a  ,   t(  ,    and  B  are  constamt  scale 
factors.   The  symbol  t  represents  the  computer  time: 
the  time  for  the  phenomena  to  occur  on  the  computer. 
The  distcince  traversed  from  the  plate  is  directly 
related  to  the  computer  time,  t  •   This  relationship  is 
presented  in  the  next  section. 

In  attempting  to  find  a  solution  to  the  boundary 
value  problem  presented  above,  it  is  assumed  that  a 
solution  to  the  boundary  value  problem  exists  for  every 
assumed  value  of  <5.   Further,  it  is  assumed  that  the 
above-mentioned  solutions  are  selectable  from  a  set  of 
discrete  solutions  to  the  initial  value  problem  of 
Equations  (47)  euid  (48)  having  the  initial  conditions 

T=0     0=0     *=0    -a.<$>oo    y=0.  (51) 

In  order  to  check  the  above  hypothesis.  Equations 
(47)  and  (48)  are  programed  for  the  auialog  computer  for 
values  of  a  =  1,  <  =   1,    and  g  =  1,  which  imply  for  this 
first  program  that 

C=T      ♦=*      ii  =    "V    .  (52) 

The  analog  momentum  equations  become 

.2     .  ., 
♦  =   ♦   -2**   -1.  (53) 

y   =   -2*$-*.  (54) 
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Figure  3  shows  the  analog  diagram  for  Equations  (53) 
and  (54) .   The  set  of  discrete  initial  conditions  on  ♦ 
is  selected  over  a  sufficiently  wide  range  to  give  a 
reasonable  topology  of  the  total  set  of  solutions  to 
this  particular  initial  value  problem. 

The  circuit  shown  In  Figure  3  is  programed  on 
the  Applied  Dynamics  Analog  Computer.   In  programing 
this  particular  computer,  the  voltage  to  the  quarter- 
square  multipliers  must  be  greater  than  10  volts  in 
order  to  keep  the  specified  computer  accuracy.   For 
this  progrcim,  constant  multipliers  accomplish  this. 
The  notation  used  in  this  diagram  is  standard  analog 
notation;  however,  the  symbols  are  defined  in  Appendix 
C. 

Figure  4  represents  the  solutions  for  4»  '  .   The 
requirement  of  satisfying  the  boundary  conditions,  Equa- 
tions (36)  and  (37),  limits  the  acceptable  range  of 
*  "  (0)  to 

0  <  *  ' ' (0)  <  1.31.  (55) 

It  turns  out  that  of  the  possible  solutions  the  values 
of  physical  interest  are  near 

(|,  '  '  (0)  =  0.5.  (56) 

Figure  4  indicates  that  for  the  range  of 
interest,  the  analog  equation  requires  rescaling  so  that 
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the  analog  output  is  at  a  greater  voltage  to  reduce  the 
error  involved.   The  precision  of  the  analog  computer 
components  used  allows  an  accuracy  of  0.1  per  cent 
based  on  a  100  volt  output. 

A  trial  and  error  procedure  indicates  that  for 
the  range  of  interest  scale  factors  of  (J  =   0.1,  K   = 
0.001,  and  /}=  0.0015  can  be  used.   As  a  result,  the 
analog  momentum  equation  takes  on  the  following  form: 

The  Initial  conditions  on  |  are  found  from 
relationship  (44) .  The  four  sets  of  initial  conditions 
used  are 

^"    =  0.5  (|)=  5.0 

=  0.4  =  4.0 

=    0.2  .  2.0   . 

The  analog  diagram  for  Equations  (57)  and  (58) 
is  shown  in  Figure  5.   The  program  results  are  shown  in 
Figures  6,  7,  and  8,  which  represent  |,  ^,    and'j/, 
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respectively.   The  results  are  as  expected.   The 
parameter,  '^^,    which  is  directly  related  to  u,  starts  at 
0,  goes  to  a  maximum  value  near  the  center  of  the  steaun 
gap,  amd  then  returns  to  zero  at  the  surface  of  the 
droplet.   The  parameter, ^,    which  is  directly  related 
to  w,  starts  at  zero  at  the  plate  and  then  reaches  its 
maximum  value  at  the  bottom  of  the  droplet.   The  param- 
eter, i,  which  is  directly  related  to  the  total  pressure, 
goes  from  a  maximum  value  of  zero  at  the  plate  to  a 
minimum  value  in  the  center  of  the  channel.   This  is 
because  the  total  pressure  head  at  the  plate  is  par- 
tially converted  into  a  velocity  head  in  the  center  of 
the  gap  between  the  plate  and  the  water  droplet,  thereby 
giving  rise  to  a  lower  pressure.   The  pressure  then 
returns  nearly  to  its  plate  value  at  the  surface  of  the 
droplet.   There  is  a  negligible  deviation  due  to  the 
small  axial  velocity  head,  which  for  this  problem  is 
negligibly  small  compared  to  the  radial  velocity  head. 
Therefore,  the  pressure  distribution  at  the  surface  of 
the  droplet  is  taken  to  be  of  the  form 

P(r,«  )  =  i  a2  -£_  (r  2  _  r^)  .       (60) 
2    ^c 

The  phenomena  discussed  above  can  be  found  in  many 


32 
textbooks  under  the  heading  of  Bernoulli's  Equation 
(15,  p.  114). 

Table  1  lists  the  Important  numerical  values  of 
the  end  points  of  Figures  6  and  7 .   The  next  section  is 
concerned  with  the  determination  of  the  velocities  and 
evaporation  rates  from  the  data  listed  in  Table  1. 

Steam  Velocities 
The  radial  and  axial  velocities  and  the  gap 
thickness  are  determined  from  the  analog  parameters 
listed  in  Table  1.   These  parameters  are  directly  re- 
lated to  the  cixial  velocity,  w,  by  Equations  (13),  (27), 
and  (41)  resulting  in 

w  =  -  2  VTV  K^.  (61) 

The  value  of  w  at  the  surface  of  the  droplet  is  given  by 

w((5)  =  -  2  VHz/C^^j.  (62) 

where  ^s   is  the  value  of  M-^  at  the  surface  of  the  drop- 
let.  The  value  of  -^x  is  tabulated  in  Table  1.   The 
parameter,  C),  is  directly  related  to  the  velocity,  u, 
by  Equations  (14),  (30),  and  (43),  resulting  In 


=  r  a^  $.  (63) 

6 
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TABLE    1 
ANALOG   CX)MPUTER   RESULTS 


$(0) 
volts 


r.=  ^-\^6      $.=  (A-VT^ 


seconds 


volts 


4(^ 


I  2 
volts 


5.0 
4.0 
3.0 
2.0 


9.9 
7.0 
5.3 
3.6 


82.0 

29.5 

13.0 

4.0 


12.25 
6.25 
3.50 
1.50 


*  The  output  of  the  ^nalog  computer  is  read  in 
volts;  however,  the  output,  cp,  is  considered  to  be 
volts  per  unit  time  when  used  in  the  equations,  in 
order  that  the  units  will  be  consistent. 
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In  a  similar  manner,  the  time  required  to  satisfy 
the  boundary  condition  of  4)  =  0,  Equation  (50),  that  is, 
the  time  at  which  the  curves  in  Figure  6  taXe  to  reach 
their  zero  values,  relates  directly  to  the  gap  thickness, 
O,  by  Equations  (26)  and  (40) ,  resulting  in 


6  =  a/^  cf  r,  .  (64) 


The  parameter,  a,  still  an  unknown,  is  now  determined 
from  the  static  equilibrium  condition.  Equation  (5) . 
Substituting  Equation  (60)  into  Equation  (5)  auid  solving 
for  a2 

2    1   Pd  n           1 
a2  =  -  V  —  I  g  ,        (65) 


I  i^o     ~   r^)r  dr 


where 


9o 


Performing  the  required  integration  and  solving  for  a, 


4/ir       V   I** 


a  =  .      ,  .^ 


o 
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Therefore,  Equations  (61),  (63),  and  (64)  taXe 
on  the  following  forms: 


u  =  r 


K  i)l±  ^  Tg  — 


o   / 


(67) 


w 


-..2K(pl±^yr.^JL\'     , 


6  = 


rt  p 


r  4 
^o 


_4_  _ft_ 


Hg. 


V 


•o  / 


(68) 


<69) 


For  a  flat  spheroid  geometry,  simplifying 
Equation  (66)  by  relating  the  radius  to  the  volume 
conveniently  results  in 


V  =  A  t  , 


(70) 


where 


A  =  TTr  2 
o 


Therefore 


4  _ 


TT^'-^ 


(71) 
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Substituting  into  Equation  (66)  results  In 


or 


P        V 


(72) 


=  vi^  rM§ 


/I 
V  / 


(73) 


Therefore   the   equations    (61),     (63),    and    (64)    take  on 
the   following  more  convenient   forms  using  the   above 
substitutions, 


u  = 


•VTrfg.4  ■^'t'' 


1 
1  V    / 


(74) 


w  = 


6  = 


-  (64rrg,)^K'C|)rM^^       ^ 


I  V  / 


,(75) 


(76) 


h  rh  „h 


^^TT^cPt)     ^     ' 


The  next  section  deals  with  the  determination 
of  the  evaporation  rate  from  the  above  velocity 
relationships . 
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Mass  Flow  Rate 
The  mass  loss  required  to  satisfy  the  condition 
of  static  equilibrium  is  calculated  by  the  relationship 


dt 


w(6) 


A  . 


(77) 


Substituting  Equations  ^70)  and  (75)^  into  the  above 
equation  results  in 


dM 
dt 


=  (64  7rgJV4)^^(A3>^(/^^^)''  ^-  OQ] 


Or  expressing 


P  = 


_1_ 

u 

Ji 

p 


(79) 
(80) 


Equation  (78)  takes  on  the  form 


^=te4rr,Acf^r'«(^)^(^) 


^  vi 


(81) 


In  this  analysis,  desiring  to  determine  the 
evaporation  rate  as  a  function  of  the  plate  temperature, 
the  temperature  dependent  parameters  of  density  emd 
viscosity  are  expressed  in  terms  of  the  film  temperature 
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of  the  vapor.   Figures  26  and  27  in  Appendix  D  indicate 
a  linear  relationship  between  the  viscosity  and  the 
specific  volume  with  the  absolute  temperature  of  the 
superheated  steam  in  the  range  of  interest. 
Therefore 

V   =      T'Tf  (82) 

A^  =  7^  T^  .  (83) 

Substituting  Equations  (82)  and  (83)  into  (81) 
results  in  the  following  form  of  the  evaporation  equa- 
tion; 

/       i 

Defining  the  material  quantity 

Q=  i647rg^  P^)   J-^  (3.281x10-2)^x453.6 

X  (3.531x10-^)-?  .  (85) 

Substituting  Equation  (85)  into  (84)  results  in 

-■or*-.    T^%^  <^^     .  (86) 

The  last  expression  gives  the  evaporation  as  a  function 
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of  the  material  properties, W,  the  gravitational  constant, 
I,  thickness  of  the  water  spheroid,  i,    the  absolute  film 
temperature,  T^,  and  the  volume  of  the  droplet  under 
consideration,  V. 

Equation  (86)  is  to  be  evaluated  for  different 
values  of  the  gap  thickness.   Substituting  Equations 
(79],  (80),  (82),  and  (83)  into  (76),  so  as  to  relate 
the  gap  thickness  to  temperature,  results  in  the  follow- 
ing relationship  for  the  gap  thickness: 

6  =  Ar^  r'^  r^  Tf^  v^  ,  (87) 

where 


"1   _  (SCflf)       (3.281x10*^1  x(. 3531x10-^)^  xl2 


.(88) 


The  pareimeter,  A  ,  is  a  function  of  material  only  as 
is  i'.   Equations  (86)  and  (87)  are  now  evaluated 
separately  to  determine  the  required  flow  rate  and  gap 
thickness  which  satisfy  the  momentum  equation,  contin- 
uity equation,  and  the  condition  of  static  equilibrium. 

What  remains  now  is  to  determine  the  evaporation 
rate  from  an  energy  consideration  as  a  function  of  the 
gap  thickness.   Therefore,  the  next  two  sections  deal 
with  the  solution  of  the  energy  equation  euid  the 
macroscopic  energy  balance. 
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Energy  Ecmation 
This  section  is  concerned  with  the  calculation 
of  the  amount  of  heat  transferred  from  the  plate  to 
the  water  droplet  by  conduction.   Previous  work  (4), 
C5)  ,  cind  (6)  used  the  relationship 

,  ,  l!^  {Tp  -  T^^^)  ,  (89) 

for  the  calculation  of  the  amount  of  heat  transferred 
through  the  gap.   However,   the  above  relationship  does 
not  consider  the  effects  of  the  stream  velocities  on 
the  conduction  heat  trcinsfer.   Therefore,  the  above 
equation  is  considered  a  first  order  approximation  to 
the  energy  equation,  Equation  (4) .   When  the  velocity 
effects  are  neglected,  the  energy  equation  takes  on  the 
form  of  the  Laplace  equation,  of  which  Equation  (89)  is 
the  solution. 

As  a  result  of  the  work  of  the  previous  two 
sections,  the  effect  of  the  stream  velocities  on  the 
rate  of  heat  transfer  by  conduction  through  the  steam 
gap  can  be  determined.   The  linear  relationship  implied 
by  Equation  (89)  is  perturbed  by  the  ejection  of  the 
saturated  steam  into  the  vapor  stream. 
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For  the  problem  under  consideration,  the 
physical  conditions  indicate  that 


—  ^  ^  —   .  (90) 

dr     Oz 


Therefore,  Equation  (4)  can  be  simplified  to 


dT     d^T 
w  —  =a^—   -  (91) 

^2     dz^ 


Substituting  relationship  (61)  into  Equation 
(91)  results  in 

, r  dT     d^T 

-  2  YH;  ;CCp— =CX —  .         192) 

^  dz     dz2 
However,  combining  Equations  (26)^  and  (40)  results  in 

Differentiating  the  above  with  respect  to  z  yields 

dT   1 


7t-  VTT   •  ^^^^ 

dz   0  V  i; 


Therefore,   the  first  and  second  derivatlvas  with  respect 
to  temperature  are  written  as  follows: 
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dz    dTdz    (f      ^  V     dT 


d^T     1   a   d^T 


dz 


2     cf'  y  .T^ 


(96) 


Substituting  the  above  two  relationships  into  Ecjuation 
(92)  results  in  the  following  form  of  the  energy  equa- 
tion: 


^  =   0  .  (97) 


d^T  ^  2VKCl^ 

dr2   ""     oc       dT 

The  parameter,  (|),  in  the  above  equation  is  a 
function  of  T,  and  is  conveniently  approximated  by  the 
following  form  (see  Figure  7) : 


.  *. 


d)  s   I?  r.  (98) 


Therefore,  the  energy  equation  becomes 

d2T       iVK  (^  \r  di 


Let 


+   -'   -  ii  Zl  =   0  .      (99) 

d72       0(  r     dT 

6 


B   =  ^Xl      ,  (100) 

a  r, 
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but  the  Prandtl  number  is  equal  to 


Pr  =  —   .  (101) 

a 


Therefore,    the  constant,    B,    takes   on  the   form 

2  Pr^c  cr(|>5 

B      =      .  (102) 

^5 


and  Equation  (97)  becomes 


d^T         dT 

—  +   bT     =   0   .  (103) 

dr  dV 


The  problem  now  is  to  integrate  the  above 
differential  equation  and  apply  the  thermal  boundary 
conditions  shown  in  Equations  (49)  cind  (50)  .   Equation 
(103)  is  made  readily  integrable  by  substituting  in 

dT 
y  =  —   •  (104) 

dT 

Substituting  Equation  (104)  into  (103)  2uid  integrating 
yields 

^  =  c^  exp  {--   bT^}  .        (105) 
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Integrating  the  above  results  in 


=  J  C4  expj 


--bT2 


dT  + 


(106) 


The  value  for  the  above  integral  is  given  in  reference 
(16,  p.  303)  as 


T   =   c 


4  1  2 


"V^erffylr])   +   C5  .   (107) 


But,  reference  (16.  p.  297)  shows  that 


erf  X  = 


00 


z 


(-1)"  x2"-^l 


Vrf  n=0   n:   (2n  +  1) 


(108) 


Expanding  Equation  (108) , 


2    .     ^3   ^5 

erf  X  =  X  -  —  +  — 

VTT  I     3    10 


(109) 


Substituting  Equation  (109)  into  (107)  results  in 


T   =    c^r 


B_,2    B'' ^4 

1  -  — T  +  — T 

6       40 


+  c 


5  • 
(110) 


Defining 


B  ^  2 


^     6  « 


(111) 
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results  in 


T    =    c^r   1  -  A 


/  r  \ 


IT. 


.^/-I 


1.1  ^  ^^5 


(112V 


The  value  of  the  parameter,  /\»  is  approximately 
0.05.   Consequently,  it  follows  that 


A^  <  <  A  . 


Since  X   is  defined  for  the  domain  0  =   T  =  T    , 

0 


(113) 


1.  u .  b. 


I   f  \ 


e  ■ 


-  1  . 


(114) 


Therefore,  from  calculus  (17,  p.  129),  it  is  known  that 

00 

an  alternating  series  of  the  form  \   (-l)'^"'""^ 

IT 


a   where 


a^a^^^j^^^O  and  where  the  limit  aj^  =  0  converges  and  the 
remainder  after  n  terms  has  a  value  between  zero  and 
the  first  term  not  taken.   Consequently,  second  order 
terms  and  higher  of  A  are  neglected,  since  the  maximum 
error  in  the  resulting  series  is  less  than  0.0025. 
Therefore,  the  temperature  is  represented  by 


T  ^  c.r  (i  -  Af^ 


+  c 


5  * 


(115) 


46 

Evaluating  c^  and  c^  from  the  boundary  conditions, 
Equations  (49)  and  (50),  yields 


T   =   Tp  -  — £ 


sat' 


u-A)    r 


i-A 


h' 


(116) 


However,  in  this  particular  problem,  the 
temperature  distribution  under  the  droplet  is  not  of 
great  interest;  rather,  the  heat  flux  at  the  droplet 
interface  is  the  important  quantity.   This  is  found 
from  the  relationship 


q=  -It  A^ 
dz 


(117) 


Using  relationship  (94)  the  above  becomes 


q  =  - 


k  A 


a 


dTl 

d? 


(118) 


Differentiating  relationship  (115)  results  in  a  tempera- 
ture gradient  of  the  form 


3T 


^^P  -  ^sat^ 

r^(i  -  A) 


-3A(f 


(119) 


Therefore,  the  heat  flux  at  the  surface  of  the  water 
droplet  is  given  by 
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k  A 


q  = 


<^\-^JR 


(T   -  T   ^) 
^  p    sat 


1  -  3A 

1-  A 


(120) 


but,  substituting  in  Equation  (64)  results  in 


k  A 


6 


<^p-^aat'-^- 


4121] 


where 


A  = 


1  -  3A 

i-A 


(122) 


The  parameter , vV,  represents  a  velocity 
correction  factor  to  the  above  equation.   Bound  up  in 
this  correction  factor  is  the  consideration  that  some 
of  the  heat  leaving  the  plate  goes  into  superheating  the 
vapor  leaving  the  surface  of  the  droplet.   The  value  of 
-/v.  is  approximately  0.95;  thus,  the  velocity  correction 
factor  represents  a  5  per  cent  correction  on  the  energy 
equation. 

Next,  the  solution  of  the  energy  equation  is 
used  in  a  macroscopic  energy  balance  to  determine  the 
allowable  evaporation  rate  as  a  function  of  ($  . 
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Macroscopic  Energy  Balance 
The  amount  of  mass  transfer  from  the  water 
droplet  is  now  calculated  explicitly  as  a  function  of 
0  by  solution  of  the  macroscopic  energy  equation.   The 
macroscopic  energy  balance  for  the  water  droplet  model 
as  shown  in  Figure  2  takes  the  form 


hfg  ^   =   q^  +  q^   .  (123) 


Here  q^  is  the  conduction  energy  flux,  q  is  the  net 
radiative  energy  flux,  h   is  the  increase  in  enthalpy 

during  the  vaporization  of  one  pound  mass  of  liquid, 

d_M 
and  is  the  amount  of  liquid  vaporized  by  the  energy 

dt 

transferred  by  conduction  and  radiation  through  the 

stecun  gap. 

The  conduction  energy  flux  is  represented  by 

Equation  (121) ,  while  the  radiative  flux  Is  given 

approximately  by  the  relationship  (18.  p.  64) , 


^r  =  °o  ^e  Fa  ^  (Tp*  -  T^^^^)  .    (i24> 
From  geometric  considerations  (19,  p.  199,  formula  6), 


F^  =   1  .  (125) 
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For  the  above  geometry,  where  F^  =  1,  the  F   factor 

a  e 

which  considers  the  departure  of  the  two  surfaces  from 
complete  blackness  is  represented  as  (18,  p.  61): 

1 
F    =     .  (126) 

1.^1.-    1 
fp    ^D 

A  water  droplet  at  212  F  acts  similar  to  a  black  body, 
since  the  emissivity  of  water  at  212  F  is  0.963  (20, 
p.  47a) .   Therefore, 


Fg   =   Cp   .  (127) 


Thus,  the  radiative  flux  is  written  as 

qr   =  ^  o  ^p  A  (Tp4  _  Tg^^4,  .       (128) 

The  above  equation  is  conveniently  rewritten  in  the  form 
q^   =   ^p  A  Ft  (T   -  Tg^^)  ,  (129) 


where 


F^   -  ,        (130) 

^Tp   -   ^sat  > 


with  the  values  of  F^  available  in  the  literature  (19, 
p.  208}  as  a  function  of  the  body  temperatures. 
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The  absorption  of  some  of  the  radiative  energy 
by  the  water  vapor  is  neglected  in  this  problem  because 
of  the  small  path  length  between  the  droplet  and  the 
plate  Cli,  p.  214)  and  (21,  p.  388). 

Thus,  the  macroscopic  energy  balance  takes  on 
the  form 

Substituting  in  Equation  (70)  and  solving  for  the 
evaporation  rate  yields 

dt     £  h^g  \  «  / 

with  all  the  temperature  dependent  properties  evaluated 
at  the  film  temperature.  Equation  (8) . 

The  overall  heat  trasfer  coefficient,  U,  between 
the  plate  and  the  water  droplet  is  defined  by 

q^  +  q^   =   U  A  (Tp  -  T^^^)  .  (133) 

Comparing  the  above  relationship  to  the  Equation  (131) 
results  in  the  following  form  of  the  overall  heat 
transfer  coefficient: 

U  =  -A+  C„  Ft   .  (134) 

6         ^ 
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Substituting  the  above  into  Equation  (132)  results  in 

dM      1 

=       U  A  (T   -  Tg^t)    •  (135) 

*^^     *^fg 

Clearly,  if  the  evaporation  rate  is  calculated  by 
theory  or  experimentally  measured,  the  overall  heat 
transfer  coefficient  is  known  directly  from  the  evapor- 
ation rates  by  the  use  of  Equation  (135). 

Following  the  earlier  procedure  of  expressing 
the  transport  parameters  in  terms  of  temperature,  the 
thermal  conductivity  is  expressed  as  a  linear  function 
of  temperature  of  the  form 

k  =  m  +  n  Tf  ,  (136) 

as  shown  in  Figure  28  of  Appendix  D. 

Therefore,  the  evaporation  rate  shown  in 
Equation  (132)  takes  on  the  form 

dM     V   /l2(m  +  nTf)A       \ 

=  +   C^Fm      X 

(■^p  -  "^sat^  ^^  ^    (1.076x10-3)  .  (137) 
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Graphical  Determination  of  Gap  Thickness 
and  Evaporation  Rate 

The  evaporation  of  water  vapor  from  a  water 

droplet  has  been  determined  in  two  ways.   First,  from 

the  standpoint  of  momentum  required  to  produce  static 

equilibrium  for  a  given  gap  thickness  (repeated  for 

convenience) , 

^  ^QT^-^tM*  ,    (86) 

dt  f     -^  6 

and  secondly,  from  a  standpoint  of  energy  transfer, 
dM       V   |12(m  +  nTf)A 


dt 


+  Cp  Ft  I  ^ 
'   ^fg^      ^ 

453.6  , 

(T   -  T^^.  )  X  (1.076x10-3) 

P    s^t   3500 

(137) 

Both  equations  shown  above  are  solved  explicitly  for 
different  values  of  the  gap  thickness  found  from  the 
relationship: 

6  =  A  "^  r~^  l-^   T^l  V^  .   (87) 

When  the  evaporation  rates  calculated  from  Equations 
(86)  and  (137)  are  equal  for  a  given  value  of  the  gap 
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thickness,  as  calculated  from  Equation  (87),  the 
governing  equations  and  boundary  conditions.   Equations 
(1)  through  (7),  are  satisfied  concurrently. 

For  example.  Figure  9  shows  a  graphical  solution 
of  the  momentum  and  energy  equations  for  the  case  of  a 
0.5  cc  water  droplet  resting  on  a  plate  at  600  F.   The 
point  of  intersection  of  the  two  equations  on  Figure  9 
represents  the  conditions  where  all  the  governing  equa- 
tions are  satisfied  concurrently.   The  values  of  the 
evaporation  rate,  which  is  directly  related  to  the 
overall  heat  transfer  coefficient  by  Equation  (135) , 
and  the  gap  thickness  are  found  directly  from  the 
ordinates  in  Figure  9.   Similar  graphs  were  constructed 
for  different  size  droplets  and  for  various  plate 
temperatures  to  determine  the  evaporation  rates  as  well 
as  the  overall  heat  transfer  coefficient  for  a  variety 
of  conditions. 

It  is  shown  in  a  later  section  that  the 
theoretical  results  are  in  excellent  agreement  with 
theory. 
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Fig.  9.   Graphical  Simultaneous  Solution  of  Momentum  and 
Energy  Equations  for  V=0.5  cc,  T  =600  F,  and  £"  =0.5. 


CHAPTER  III 

FLOW  DISTRIBUTION 
The  velocity  distribution  and  the  Reynolds 
number  under  the  droplets  are  considered,  since  the 
earlier  assumption  of  laminar  flow  under  the  droplet 
is  now  justified.   In  addition,  it  is  important  to 
determine  the  magnitude  of  the  radial  velocity  at  the 
edge  of  the  droplet,  since  the  droplet  tends  to  move 
slowly  about  when  resting  on  a  flat  plate. 

The  velocities,  u  amd  w,  are  related  to  the 
computer  variables,  *  and  *  ,  by  Equations  (61)  and 
(63).   For  a  given  set  of  physical  conditions,  such 
as  volume  of  the  droplet  and  temperature  of  the  heat- 
ing plate,  the  velocity  distributions  vary  directly 
as  *  and  * .   Consequently,  the  curves  shown  in  Figures 
6  and  7  represent  the  forms  of  the  radial  and  axial 
velocity  distributions.   Previously,  the  exact  solu- 
tion for  the  velocity  distribution  could  not  be 
determined  because  the  gap  thickness,  6  ,  was  an  unknown. 
However,  using  the  results  of  the  previous  section, 
the  velocity  distribution  can  be  evaluated  since  the 
gap  thickness  is  now  a  known  quantity. 
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The  curves  shown  in  Figures  6  and  7  are 
generated  from  known  values  of  (I)(0)  rather  than  specify- 
ing 0.  However,  0  is  related  to  0(0)  by  relationship 
(87)  and  the  curve  shown  in  Figure  10  is  constiructed 
from  this  relationship.   Thus,  for  a  known  0,  Cp(0)  can 
be  determined  directly  from  Figure  10.   However,  a  close 
inspection  of  Figure  6  indicates  that  for  all  practical 
purposes  the  curves  representing  (D  are  parabolas  which 
can  be  fitted  by  the  equation 

2r      i2  /    (i) 

—  -  1    =   -  —^ 1  K  (138) 


-^1  2 


or  using  relationship  (64) 
(2z     ■' 


^ 


—  -  1    =   -  — ^ 1    .  (139) 

0       I  \^n^ 


The  above  equation  can  be  rearranged  to  the  form 


(140) 


The  values 


°4{^) 


in  the  above  ec[uation  are 


plotted  in  Figure  11  as  a  function  of  6.   This  figure 
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Fig.    10.      0  as   a  Function  of  (^(0)    for  v=0.5  cc, 


T   =600  F,    and  C  =0.5. 
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Is  constructed  from  the  analog  results  presented  in 
Figure  6  and  the  information  presented  in  Figure  10. 
The  vertex  of  the  parabola  is  taken  as  the  anchor 
point  between  the  analog  data  and  Equation  (140) . 

The  parabolic  relationship  very  accurately 
represents  the  curves  in  Figure  6  as  seen  in  Table  2, 
which  presents  a  comparison  of  the  C^  calculated  from 
Equation  (140)  and  from  the  actual  data  shown  in 
Figure  6,  for  the  specific  case  of  <^(0)  =  4.0.   Thus, 
the  radial  velocity  distribution  across  the  steam  gap 
can  be  considered  to  be  parabolic  in  shape.   Conse- 
quently, for  parabolic  flow  the  average  velocity  is 
two-thirds  the  maximum  velocity  (22,  p.  62K   There- 
fore the  average  radial  velocity,  as  defined  in  Equa- 
tion (9) ,  takes  on  the  form 

u  =  2/3  r  a  -$>[—)   •  (141) 

The  maximum  radial  velocity  occurs  at  the  edge  of  the 
droplet  where  r  =  r^;  the  maximum  average  radial 
velocity  is  expressed  as 

%=  2/3  r^a-C^^I   .  (142) 
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TABLE    2 


PARABOLIC    FIT    OF  ^ANALOG    RESULTS    FOR  Cp(0)    =4.0 


^  $cal=  ^Fi 


g.    6 


0 

0.0 

0.0 

0.0 

6/8 

0.875 

2.69 

2.7 

6/4 

1.75 

4.59 

4.6 

36/8 

2.625 

5.75 

5.75 

6/2 

3.5 

6.125 

6.12 

56/8 

4.375 

5.75 

5.7 

36/4 

5.25 

4.59 

4.55 

76/8 

6.125 

2.69 

2.7 

6 

7.0 

0.0 

0.0 
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The  meiximum  Reynolds  number  beneath  the  droplet  is 
calculated  from  the  flow  beneath  the  droplet  by  the 
relationship 

Re  =  .  (143) 


where 


^       .    flow  cross  section       ,, ..x 
Dg  =   4  X  ,     (144) 

wetted  perimeter 


=  4  X 


=  26. 


2'^r^6 


2  X  2TrQ 


Thus 


v^r^ 


r 


R 


=  a/3  L£ll_Vr  :56U-  .(145) 


Consider,  for  example,  the  previous  problem  of 
a  0.5  cc  droplet  on  a  600  F  flat  plate.   Figure  9 
indicates  that  the  gap  thickness  is  0.00475  inches. 
From  Figure  11,  the  value  of  cp 


is  equal  to  13.25 


2  / 
Using  these  values  in  Equations  (142)  and  (145) , results 

in  a  Reynolds  number  of  10.6  and  an   average  radial 
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velocity  at  the  edge  of  the  droplet  of  5.25  ft/sec. 
Thus,  the  flow  is  well  within  the  laminar  range  smd  the 
slight  motion  of  the  water  droplet  on  the  heating  plate 
is  small  compared  to  the  average  steam  velocity  leaving 
the  gap  beneath  the  droplet. 

An  increase  in  the  volume  of  the  droplet 
increases  the  exit  radial  steam  velocity  and  the  Rey- 
nolds number,  since  both  ^  and  <5  increase  with 
increasing  volume.   Also,  combining  Equations  (74), 
(79),  and  (32)  indicates  that  an  increase  in  temperature 
of  the  heating  plate  increases  the  exit  radial  steam 
velocity  and  the  Reynolds  number  in  proportion  to  (Tf)^. 
However,  for  the  temperature  range  and  volume  range 
investigated  in  this  paper,  the  basic  conclusion  that 
the  flow  is  laminar  is  not  affected  by  the  volxime  and 
temperature  changes  considered. 


CHAPTER  IV 

EXPERIMENTAL  PROCEDURES 

In  the  theoretical  analysis,  the  mass  evaporation 
rate  for  a  droplet  on  a  flat  plate  is  shown  to  be  a 
function  of  the  plate  temperature  and  the  volume  of  the 
droplet.   Hence,  an  experimental  verification  of  the 
theory  requires  that  the  evaporation  rate  be  measured 
for  different  plate  temperatures  and  droplet  volumes. 

The  evaporation  rate  is  determined  experimentally 
from  measurements  taken  on  the  total  vaporization  time. 
The  total  vaporization  time,  that  time  required  for  the 
entire  volume  of  liquid  which  is  placed  on  a  heating 
surface  to  vaporize  completely,  is  measured  as  a  func- 
tion of  droplet  size  for  various  plate  temperatures  and 
surface  conditions.   The  experimental  data  are  listed 
in  Table  6  of  Appendix  E,  while  the  plots  of  the  data 
are  shown  in  Figures  12,  13,  and  14.   The  slopes  of 
these  curves,  rate  of  change  of  volume  with  respect  to 
time,  represent  the  evaporation  rate  of  the  droplet. 
The  determination  of  the  slope  of  a  curve  V  =  f(t),  when 
a  table  of  distinct  sets  of  values  (Vj^/t  )  are  known,  is 
considered  in  detail  in  the  next  chapter. 
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Fig.  12.   Total  Vaporization  Time  for  Water  Droplets 
on  a  Flat  Plate  as  a  Function  of  Their  Initial 
Volume  for  Various  Surface  Conditions  at  a 
Plate  Temperature  of  Approximately  600  F. 
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of  the  Heating  Surface  which  Had  a  1°  Apex  Angle. 
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Fig.  14.   Comparison  of  the  Total  Vaporization  Time 

for  Water  Droplets  on  a  Flat  Plate  and  a  1° 

Conical  Surface  at  Approximately  600  F. 
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The  test  sections  used  for  vaporizing  the  water 
droplets  are  shown  in  Figures  15  and  16.   Both  a  304 
stainless  steel  and  a  graphite  test  plate,  as  shown  in 
Figure  15,  were  fabricated  to  allow  a  wide  variation 
in  surface  conditions.   The  stainless  steel  is  a  hard 
metal  impervious  to  the  liquid,  while  the  graphite 
exhibits  many  small  cracks  across  its  surface.   The 
heating  surfaces  of  the  plates  were  machined  to  a 
satin  finish.   A  satin  finish  is  equivalent  to  a 
surface  finish  of  approximately  125  microinches,  rms . 
A  depth  of  cut  of  0.001  inches  with  a  cutting  speed 
of  0.0014  inches  per  revolution  was  used.   The 
machined  surfaces  were  flat  to  0.0005  inches,  as  veri- 
fied by  use  of  a  dial  indicator  gauge  mounted  on  the 
lathe  carriage  which  had  performed  the  finishing  cut. 
In  addition  to  the  data  taken  with  the  above  surfaces, 
some  data  on  a  flat  polished  brass  surface  were  avail- 
able from  reference  C5)  . 

The  data  shown  in  Figure  12  indicate  quite 
plainly  that  the  surface  condition  has  no  noticeable 
effect  on  the  vaporization  time.   The  volume  range 
below  1.5  cc  in  Figure  12  represents  the  small  sphe- 
roidal and  flat  spheroidal  region,  while  for  initial 
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Fig.  15.   Schematic  Cross  Section  of  304  ss 
and  Graphite  Test  Plate. 
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Fig.  16.   Schematic  Cross  Section  of  Test  Plate  with 
a  1  Conical  Heating  Surface. 
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volumes  greater  than   1.5  cc  the  droplet  Is  in  the  bubbly 
spheroidal  region,  as  depicted  in  Figure  1. 

The  surface  condition  does  not  noticeably  affect 
the  vaporization  time  in  either  the  flat  spheroid  or 
the  bubbly  spheroidal  region,  because  the  droplet  is 
resting  on  its  own  vapor  which  prevents  an  interaction 
of  the  surface  with  the  droplet.   This  agrees  with  the 
previous  experimental  observations  (23,  p.  191)  that 
the  friction  factor  is  independent  of  the  surface 
roughness  under  conditions  of  laminar  flow. 

A  slight  problem  with  the  flat  heating  surface 
results  from  the  movement  of  the  water  droplet  against 
the  barrier  wall  during  the  vaporization  process.   How- 
ever, the  interaction  between  the  wall  with  the  water 
droplets  in  the  volume  range  of  interest,  1  cc  or  less, 
is  negligibly  small,  since  only  a  very  small  fraction 
of  the  droplet 's  periphery  touches  the  barrier  wall. 
Nevertheless,  to  eliminate  the  effect  of  contact  with 
the  barrier  wall  on  the  experimental  evaporation  rate, 
the  experimental  data  to  be  used  in  comparison  with 
the  theoretical  results  were  taken  on  a  test  section 
with  a  1  degree  apex  angle,  as  shown  in  Figure  16. 
Figure  14  presents  a  comparison  of  the  total  vaporiza- 
tion times  as  measured  on  a  flat  surface  and  a  conical 
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surface  with  a  1  degree  apex  angle.   As  seen  in  this 
figure,  there  is  no  noticeable  difference  in  the 
vaporization  times,  thereby  confirming  the  earlier 
observations  that  the  side  wall  interaction  is 
negligibly  small. 

The  test  sections  were  mounted  on  the  base 
plate  shown  in  Figure  17.   Pyrex  brand  wool  amd  glass 
wool  insulation  were  packed  around  the  heating  plate 
to  reduce  the  heat  loss  from  the  ends  of  the  test 
plate,  thereby  giving  a  more  uniform  temperature  dis- 
tribution across  the  test  plate.   A  three  inch  aluminum 
fence  was  built  around  the  test  section  to  reduce  the 
convection  currents  that  would  tend  to  flow  over  the 
surface  because  of  the  free  convective  heat  loss  from 
the  test  plate.   In  addition,  the  fence  more  nearly 
produces  the  condition  in  which  the  water  droplet  is 
surrounded  by  saturated  vapor,  thereby  reducing  any 
mass  transfer  from  the  top  of  the  droplet. 

The  test  plate  and  base  plate  were  mounted  on 
a  1200  watt,  220  volt  electrical  heating  unit.   The 
amount  of  current  to  the  heating  unit  was  controlled  by 
a  variac  as  shown  in  Figure  18. 

Pour  20  gauge  Qiromel-Alumel  thermocouples  were 
embedded  beneath  the  test  section  at  positions  indicated 
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220  Volts 


Four  Chrome 1-Alumel  Thermocouples 


Reference  Junction 
in  Ice  Bath 


Potentiometer  -  Honeywell  Rubicon 
(Manually  Balanced) 


Fig.  18.   Schematic  Diagram  of  Experimental 

Apparatus . 
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in  Figures  15  and  16.   The  surface  temperatures  were 
found  by  linearly  extrapolating  the  upper  and  lower 
thermocouple  reading  at  the  center  point  of  the  heat- 
ing plate  to  the  surface  of  the  plate.   The  linear 
correction  applied  to  the  center  thermocouple  was  also 
applied  to  thermocouples  near  the  edge  of  the  plate. 
The  thermocouples  were  fed  through  a  selector  switch 
to  a  Honeywell  manual  potentiometer.   A  32  F  ice 
reference  junction  was  used. 

Distilled  water  at  its  saturation  temperature 
was  placed  onto  the  surface  by  means  of  calibrated 
pipettes .   A  1  ml  and  a  10  ml  pipette  were  used  in  the 
experiment.   The  1  ml  and   10  ml  pipette  were  read  to 
an  accuracy  of  0.005  ml,  and  0.05  ml  respectively. 
Although  the  pipettes  are  calibrated  for  a  liquid  at 
20  C,  heating  the  pipettes  to  the  saturation  tempera- 
ture of  water  (100  C)  does  not  affect  the  accuracy  of 
the  volume  measurements  due  to  the  relatively  small 
amount  of  volumetric  thermal  expansion  involved.   The 
change  of  volume  due  to  an  increase  in  temperature  can 
be  estimated  from  the  relationship  VB^,   T.   The  volu- 
metric expansion  estimated  by  the  previous  expression 
affects  the  results  only  if  it  were  possible  to  measure 
the  volume  to  four  significant  figures. 
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The  times  in  which  the  distilled  water  is  ejected 
from  the  pipette  (shown  in  Figure  19)  are  short  and  do 
not  significantly  affect  the  vaporization  curves.   How- 
ever, because  the  ejection  time  is  short,  the  water 
leaves  the  pipette  in  a  fine  jet,   When  this  jet  was 
allowed  to  impinge  directly  onto  the  surface  of  the 
heating  plate  at  600  F  or  onto  the  top  of  a  droplet 
resting  on  a  600  F  plate,  the  cooling  effect  of  the  jet 
at  the  point  of  surface  contact  (see  Figure  20,  a  and 
b)  initiates  nucleate  boiling  which  evaporates  a  con- 
siderable amount  of  liquid  in  a  very  short  time. 
Consequently,  the  jet  was  prevented  from  impinging  di- 
rectly onto  the  heating  surface  by  directing  the  jet 
against  a  non-wettable  surface  in  the  manner  shown  in 
Figure  20  c.   When  the  heating  plate  temperature  was 
set  at  1000  F,  the  cooling  effect  of  the  jet  did  not 
initiate  nucleate  boiling;  consequently,  it  was  not 
necessary  to  use  the  technique  shown  in  Figure  20  c 
at  this  higher  temperature. 

In  collecting  the  data,  the  variac  was  first  set 
at  a  desired  value  and  the  equipment  was  allowed  to  warm 
up  slowly  to  a  steady  state  value.   Normally,  this  re- 
quired 2  to  3  hours  depending  on  the  required  surface 
temperature.   When  a  run  was  made,  the  thermocouple  emfs 
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Fig.  19.   Ejection  Time  of  Water  Droplet  from 
Pipette  to  the  Hot  Plate  Surface  as  a  Function 
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were  usually  recorded  from  the  Rubicon  before  the 
droplet  was  placed  on  the  plate,  some  time  during  the 
vaporization  process,  and  immediately  following  the 
vaporization.   Because  of  the  relatively  large  amount 
of  time  required  for  vaporization,  there  was  ciraple 
time  to  record  all  measurements  by  hand. 

The  results  of  the  experimental  measurements 
are  discussed  next. 


CHAPTER  V 

EVAPORATION  RATES 
Theoretical 

The  theoretical  determination  of  the  evaporation 
rates  for  various  plate  temperatures  and  volumes 
requires  the  solution  of  Equations  (86),  <87) ,  and  (137) 
along  with  the  construction  of  graphs  similar  to  that 
shown  in  Figure  9.   Equations  (86),  (87),  and  (137)  were 
programed  for  solution  on  the  IBM  709  digital  computer 
and  the  compiled  results  were  used  in  the  graphical 
solutions.   The  material  parameters  used  in  the  analysis 
are  tabulated  in  Appendixes  D  and  F. 

The  theoretical  mass  evaporation  rates  are  shown 
in  Figure  21  as  a  function  of  droplet  volume  and  plate 
temperature.   This  figure  contains  the  locus  of  the 
graphical  solution  points  (see  Figure  9)  for  various 
droplet  volumes  and  plate  temperatures .   As  seen  in 
Figure  21,  the  evaporation  rate  increases  with  increas- 
ing volume  of  the  droplet  and  with  increasing  plate 
temperature.   The  increase  in  the  evaporation  rate  with 
increasing  droplet  volume  is  due  primarily  to  the 
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Fig.  21.   Theoretical  Mass  Evaporation  Rate  of  a  Water 

Droplet  as  a  Function  of  Volume  for  a  Plate  Enunisivity 

of  0.5  and  Plate  Temperatures  of  600  F  and  1000  F. 
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increase  in  heat  transfer  area  associated  with  the 
increase  in  the  droplet  volume,  while  the  increase  in 
the  evaporation  rate  due  to  an  increase  in  temperature 
is  due  primarily  to  the  increase  in  the  thermal  con- 
ductivity of  the  steam.   The  temperature  dependence  of 
the  viscosity  and  specific  volume  has  a  relatively 
slight  effect,  as  seen  in  Equation  (34)  .   Here  it  is 
observed  that  the  temperature  effect  is  dampened  by 
the  one-quarter  power  on  the  absolute  film  temperature. 

The  calculated  gap  thickness,  as  shown  in 
Figure  22,  is  relatively  insensitive  to  volume  changes, 
but  is  affected  by  increased  plate  temperatures . 

Experimental 
In  attempting  to  compare  experiment  to  theory, 
it  is  necessary  to  determine  the  slope  of  a  curve 
V  =  f (t)  prescribed  by  a  set  of  tabulated  values  (V^, 
t^) .   Pictorially,  a  tamgent  line  is  constructed  to  a 
graphically  fitted  curve.   However,  Lipka  (24,  p.  234) 
points  out  that  exact  or  even  approximate  construction 
of  a  tangent  line  to  a  curve  is  difficult  and  inaccu- 
rate.  Reilly  (25)  suggests  finding  the  slope  of 
tabulated  data  by  differentiating  a  polynomial  fit  of 
the  data. 
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Fig.  22.   Gap  Thickness  of  the  Water  Droplet  as  a  Function 
of  Volume  for  a  Plate  Emissivity  of  0.5  and  Plate  Tempera- 
tures of  600  F  and  1000  F. 
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The  total  vaporization  time,  listed  in  Appendix 
E,  is  used  to  determine  a  third  order  polynomial  fit. 
The  fit  is  performed  on  the  IBM  709  digital  computer 
using  the  UF-NILLS  code  which  is  currently  in  use  by 
the  Nuclear  Engineering  Department  of  the  University 
of  Florida.   The  polynomial  is  of  the  form 

V  =  P(l)t  +  P(2)t2  +  P(3)t3,  (146) 

where  t  is  the  time  required  to  completely  vaporize 
a  droplet  of  initial  volume,  V-   The  coefficients  in 
Equation  (146)  are  listed  in  Table  3  as  a  function  of 
the  plate  temperature.   They  were  determined  by  minimizing 
the  weighted  squares  of  the  residuals.   Weights  of  V~ 
%*ere  used. 

TABLE  3 
POLYNOMIAL  COEFFICIENTS 

Temperature 

of  Plate  (F)     P(l)  P(2)  P(3) 

608       -7.2266-10~5   +4.8950-10"^   +4.4080-10"^ 
1014       +4.0022-10"5   +1.0295-10"^   +5.4745-10"® 

The  curves  shown  in  Figure  13  are  drawn  from 
Equation  (146)  using  the  coefficients  listed  in  Table 
3.   The  mass  evaporation  rates  are  now  determined 
directly  from  the  polynomial  equation  (146)  by 
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differentiating  it  with  respect  to  time  and  by 
multiplying  it  by  the  density  of  the  droplet.   Thus, 

—     =     P     [p(l)  +  2Pt2)t  +  3P(3)t2  1.  (147) 
dt        ^  ^ 

The  evaporation  rates  as  calculated  from  the 
above  equation  are  plotted  in  Figure  23  as  a  function 
of  initial  droplet  volume  and  plate  temperature  along 
with  the  theoretical  evaporation  rates . 

Comparison  of  Experiment  to  Theory 
The  theoretical  and  experimental  results  are 
shown  jointly  in  Figure  23.   The  emissivities  chosen 
in  the  theoretical  calculations  are  based  on  data 
tabulated  in  reference  (20,  p.  475) .   As  seen  in 
Figure  23,  excellent  agreement  exists  throughout  the 
volume  and  temperature  range  considered.   The  devia- 
tion of  theory  and  experiment  is  less  than  5  per  cent 
in  the  volume  range  of  0.5  to  1  cc,  while  approxi- 
mately 20  per  cent  at  a  droplet  volume  of  0.05  cc. 

The  deviation  seen  at  the  lower  droplet 
volumes  is  probably  a  result  of  the  increased  devia- 
tion of  the  flat  spheroid  model  from  the  actual 
physical  situation.   The  droplet  has  a  greater  tendency 
towards  a  spherical  shape  at  these  lower  volumes. 
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Fig.  23.   Theoretical  and  Experimental  Mass  Evaporation 
Rates  of  Water  Droplets  as  a  Function  of  Droplet  Volume, 
Plate  Temperatures  and  Plate  Emissivity. 
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During  the  vaporization  process  the  droplet 
tends  to  vibrate,  particularly  at  the  higher  plate 
temperature.   The  natural  modes  of  vibration  of  the 
water  droplets  are  similar  to  those  exhibited  by  a 
ring  in  transverse  vibration.   However,  in  the  case 
of  a  liquid  water  droplet  vibrating,  the  droplet  has 
a  tendency  to  oscillate  at  one  of  its  harmonics  with 
the  other  harmonics  being  suppressed.   The  second, 
third,  and  fourth  harmonic  were  observed  in  the  exper- 
iments.  On  the  basis  of  the  agreement  of  experiment 
and  theory,  it  appears  that  the  oscillations  do  not 
significantly  affect  the  rate  of  evaporation. 


CHAPTER  VI 

OVERALL  HEAT  TRANSFER  COEFFICIENTS 
The  overall  heat  transfer  coefficient,  as 
defined  by  Equation  (134) ,  is  now  determined  from  the 
theoretical  evaporation  rate.   Solving  Equation  (135) 
for  U, 

dM 
h^  — 
fg  dt 
U  =  .  (148) 

^  (Tp  -  -^sat^^ 

where  the  area  is  expressed  in  terms  of  the  droplet 
volume  and  thickness,  Equation  (70) .   The  above  equa- 
tion is  evaluated  for  various  droplet  volumes  and   plate 
temperatures  of  600  F  and  1000  F  using  the  theoreti- 
cal mass  evaporation  rates  shown  in  Figure  21.   The 
results  are  shown  in  Figure  24.   In  the  particular 
examples  shown  in  Figure  24,  it  is  interesting  to 
note  that  there  is  a  decrease  in  the  overall  heat 
transfer  coefficient  for  an  increase  in  temperature. 
This  is  a  result  of  the  manner  in  which  U  is  defined. 
Although  the  mass  evaporation  rate  is  higher  at  1000  F 

87 


88 


fa 

CM 

m 
u 


uu 

1           ' 

1            1            I            1            ,      ■     -, 

I 

- 

- 

80 

rl 

- 

"  V\ 

- 

60 

\     \ 



i                ^^^^^^"^ 

^-..^^^                     Tp=600    F 

L 

^~^^^^==ir:^^ — 

- 

40 

-Tp=1000   F 

- 

20 

- 

- 

n 

1            1 

1             1             1             1             1             1 

1 

0.2 


0.4 


0.5 


0.8 


1.0 


V  (cc) 


Fig.  24.   Theoretical  Heat  Transfer  Coefficient  of  a  iVater 
Droplet  as  a  Function  of  Volume  for  Plate  Temperatures 
of  600  F  and  1000  F  for  a  Plate  Emissivity  of  0.5. 
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than  at  600  F,  the  temperature  difference  in  Equation 
<148)  is  also  higher  resulting  in  a  net  decrease  in 
the  overall  heat  transfer  coefficient. 

The  overall  heat  transfer  coefficient  increases 
at  lower  volumes  as  a  direct  result  of  the  lower  gap 
thickjiess  associated  with  the  smaller  droplet  volumes. 
However,  there  is  a  net  decrease  in  the  mass  evapora- 
tion rate  at  these  lower  volumes,  as  shown  in  Figure 
21,  because  of  the  smaller  heat  transfer  areas  associ- 
ated with  these  low  volumes . 

The  amount  of  thermal  radiation  is  calculated 
to  be  2  (BTU/hr  ft2}  at  600  F  and  4.4  at  1000  F. 
Consequently,  radiation  heat  transfer  at  600  F  repre- 
sents less  thcui  5  per  cent  of  the  overall  heat 
transfer  coefficient;  while  at  1000  F  it  represents 
less  than  10  per  cent.   As  the  droplet  volume  decreases, 
the  percentage  of  radiative  heat  transfer  decreases 
still  further  because  of  the  relative  increase  of 
conduction  heat  transfer  which  results  from  the  small 
gap  thickness  at  the  lower  droplet  volumes.   Therefore, 
in  the  temperature  range  under  consideration  in  this 
dissertation,  the  results  are  relatively  insensitive 
to  the  radiative  parameters  of  the  system. 
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The  magnitude  of  the  heat  transfer  coefficient 
is  one  to  two  orders  of  magnitude  below  the  heat  trans- 
fer coefficients  associated  with  nucleate  boiling 
which  occurs  when  the  temperature  of  the  heated 
surface  is  slightly  higher  than  the  saturation  temper- 
ature of  the  liquid.   The  low  magnitude  of  the 
coefficient  accounts  for  the  relatively  long  lifetime 
of  the  droplet  on  the  hot  plate. 


CHAPTER  VII 

GRAVITATIONAL  EFFECTS 
With  man  soon  to  be  traveling  to  the  moon  in 
the  next  decade,  the  effects  of  reduced  gravity  on 
the  mass  evaporation  rate  were  investigated  for  the 
particular  case  of  evaporation  in  em  enclosure 
located  on  the  moon's  surface.   The  environmental 
conditions  inside  the  enclosure  are  assumed  similar 
to  atmospheric  conditions  on  earth. 

The  effect  of  the  gravitational  force  field 
on  the  mass  evaporation  rate  is  found  by  solution  of 
Equations  (86),  (87),  and  (137)  simultaneously  as 
before.   In  the  problem  under  consideration  here,  the 
value  of  I  for  the  moon's  surface  is  0.16.   However, 
one  complication  is  introduced,  in  that  ^  is  a  func- 
tion of  the  gravitational  potential. 

Den  Hartog  (26.  p.  86)  points  out  that  the 
shape  of  a  droplet  results  as  a  compromise  between  the 
forces  due  to  capillary  tension  in  the  surface  of  the 
drop  and  the  forces  of  gravity  on  the  liquid.   The 
capillary  tension  tends  to  maXe  a  drop  purely  spherical 
as  in  the  case  of  a  freely  falling  droplet,  while  the 
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action  of  gravity  tends  to  flatten  the  droplet  as  in 
the  case  of  oil  drops  being  placed  on  water.   The 
problem  now  is  to  calculate  the  effect  of  a  decrease 
in  the  gravitational  field  on  I . 

The  physical  situation  of  a  drop  resting  on  a 
vapor  film  is  analogous  to  the  drop-shaped  storage 
tank  in  which  only  membrane  stresses  are  assumed  to 
exist.   In  the  case  of  the  droplet,  a  thin  surface 
film  of  uniform  tension  is  formed  which  envelopes  the 
liquid  and  prevents  it  from  spreading  over  the  sup- 
porting surface.   Timoshenko  and  Woinowsky-Krieger 
(27,  p.  444  and  445)  list  the  differential  governing 
equations  required  for  the  solution  of  the  droplet 
shape.   The  solution  to  this  equation  is  not  given 
explicitly;  rather,  a  numerical  integration  is  required. 
The  differential  equations,  method  of  solution,  and 
values  of  £  on  both  the  moon's  and  earth's  surface 
are  shown  in  Appendix  F. 

Using  the  calculated  values  of  £  and  I  =0.16, 
the  mass  evaporation  rate  is  calculated  in  the  manner 
described  previously.   Figure  25  shows  the  mass  evapor- 
ation rate  on  the  moon  for  a  plate  temperature  of  600  F 
and  an  eraissivity  of  0.5.   Also,  the  mass  evaporation 
rate  on  the  earth  for  an  identical  plate  temperature 
and  emissivity  are  shown  for  comparison  purposes. 
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Fig.  25.   Theoretical  Evaporation  Rates  of  a  Water  Droplet 
in  Both  the  Earth's  and  Moon's  Gravitational  Fields  for  a 
Plate  Temperature  of  600  F  and  an  Emissivity  of  0.5. 
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As  seen  in  this  figure,  the  mass  evaporation 
rate  is  approximately  half  of  that  which  exists  on 
the  earth  for  a  given  droplet  volume.   This  effect 
is  due  to  the  decrease  in  the  heat  transfer  area 
because  of  the  smaller  gravitational  field.   The 
droplet  tends  to  remain  spherical  rather  than  pancake- 
shaped  as  on  the  earth. 

It  is  hoped  that  experimental  verification 
of  the  results  shown  in  Figure  25  will  be  possible  in 
the  very  near  future. 


CHAPTER  VIII 

CONCLUSIONS 

The  mass  evaporation  rate  for  a  water  droplet 
supported  by  its  own  sup>erheated  vapor  on  a  flat 
plate  is  determined  by  the  simultaneous  solution  of 
the  momentum,  continuity,  and  energy  equations,  and 
by  experimental  measurements.   The  rate  of  vaporiza- 
tion of  a  water  droplet  is  shown  to  depend  on  the 
volume  of  the  droplet,  temperature  of  the  plate,  and 
the  gravitational  field.   The  theoretical  and 
experimental  mass  evaporation  rates  are  found  to 
agree  within  20  per  cent  over  a  volume  range  of  0.5 
cc  to  1.0  cc  and  over  a  temperature  range  of  600  F 
to  1000  F.   In  this  parameter  range,  the  mass  evapora- 
tion rate  varies  from  0.001  to  0.01  (g/sec) .   The  mass 
evaporation  rate  is  found  to  increase  for  increasee 
in  either  the  volume  of  the  water  droplet  or  the 
plate  temperature. 

Theoretically,  the  steam  flow  beneath  the 
droplet  is  shown  to  be  laminar  with  a  radial  velocity 
distribution  which  is  parabolic  in  shape,  and  the 
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steam  gap  beneath  the  droplet  is  shown  to  range  between 
0.003  and  0.008  inches.   The  gap  thickness  increases 
for  increased  volume  and  temperature,  as  did  the 
evaporation  rate.   The  shape  of  the  droplet  has  been 
calculated  by  solution  of  the  equilibrium  equations 
based  on  membrane  stress  analysis.   The  theoretical 
results  are  in  good  agreement  with  experiments. 

The  effect  of  the  gravitational  potential  on 
the  mass  evaporation  rate  is  considred  in  the  theoreti- 
cal development.   A  reduction  in  the  gravitational 
potential  from  1  (earth)  to  0.16  (moon)  reduces  the 
mass  evaporation  rate  by  approximately  half. 

A  velocity  correction  factor  to  be  used  in 
conjunction  with  Fourier's  steady  state  one-dimensional 
heat  conduction  equation  is  derived.   This  parameter 
takes  into  account  the  effect  of  the  flow  distribution 
on  the  energy  transferred  from  the  plate  to    the  drop- 
let.  In  the  temperature  and  volume  range  under 
consideration  here,  this  factor  has  a  5  per  cent  effect 
on  the  results . 

A  relationship  between  the  mass  evaporation 
rate  and  the  overall  heat  transfer  coefficient  is 
presented.   The  overall  heat  transfer  coefficient 
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ranges  between  70  (BTU/hr  ft2)  for  0.05  cc  droplets 
and  40  (BTU/hr  ft2)  for  1  cc  droplets  in  the  tempera- 
ture range  considered. 

Finally,  it  is  concluded  that  the  analytical 
model  developed  in  this  dissertation  can  be  used  to 
predict  the  evaporation  rate  and  overall  heat  transfer 
coefficients  for  water  droplets  on  a  flat  plate  in 
the  film  boiling  regime  with  a  reasonable  degree  of 
accuracy  over  the  range  of  parameters  investigated. 
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APPENDIX  A 


REACTIVE  FORCE 

The  ejection  of  mass  from  the  water  droplet 
into  the  steam  gap  below  the  droplet  produces  a 
force  on  the  droplet  itself,  much  the  seime  as  that 
produced  by  a  rocket  engine.   This  force  is  not  con- 
sidered in  the  analysis,  since  it  is  negligibly  small, 
as  is  now  shown. 

Starting  from  Newton's  law, 

1  d 
P„   = (momentum).        (149) 

""    gc  dt 

For  a  stationary  droplet  operating  under  pseudo-steady 
state  conditions,  the  above  equation  takes  on  the  form: 

w(6)   dM 
Fj^  =  —   .  (150) 

gc   <it 

The  axial  velocity,  w(6),  is  related  to  the 
evaporation  rate  by  the  continuity  equation,  which  is 
of  the  form: 


dM 

—  =  A  p  wC5  )  .  (151) 

dt 
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Therefore,  substituting  the  above  equation  back  into 
Equation  (151) ,  results  in 

2 


1      dM 

A  p  g    Idt  / 


(152) 


But,  substituting  in  Equation  (70)  into  the  above 
relationship  results  in 

2 

0      /dM 
F 


^     V  P  g 


c 


(153) 

Idt  ■ 


dM 
The  values  of  —  are  found  in  Figures  21  and  23 

dt 

while  the  other  parameters  are  found  in  Appendixes  D  and  F 

Substituting  the  numerical  values  for  the  0.5  cc 

droplet  example  problem  of  Chapter  II  into  the  above 

equation  results  in  a  calculated  reactive  force  of 

4.6  X  10"   pounds-force.   This  force  is  clearly 

negligible  when  compared  to  the  weight  of  the  water 

droplet. 


APPENDIX  B 


SOLUTION  OF  THE  MOMENTUM  EQUATION  FOR  P(r) 
In  attempting  to  solve  a  physical  problem,  it 
is  customary  to  first  use  the  simplest  analytical 
model  and  then  work  up  to  more  sophisticated  models 
in  successive  steps.   However,  if  the  pressure  is 
assumed  to  be  solely  a  function  of  the  radial  position, 
only  the  trivial  solution  for  the  velocity  distribution 
is  possible. 

Assuming  P(r),  Equation  (22)  takes  on  the  form 

2  f  f   =  -If"       .  (154, 
Differentiating  the  above  equation 

2  f  f"  +  2  f'2  =  .^f...  .  ^^35j 

Plugging  Equation  (155)  into  Equation  (21)  results  in 

3  f'^  -  a2  .  (156J 
Or,  rewriting, 

f   =   constant  .        (157) 

The  general  solution  to  the  above  ordinary  differential 
equation  is  of  course  of  the  form: 
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f  =  Cg  z  +  c^   .  (158) 

The  only  solution  to  the  above  equation  that  can 
satisfy  the  boundary  conditions,  Equations  (23)  and 
(24) ,  is  the  trivial  solution 

f  =   0  .  (159) 

It  is  really  not  surprising  that  the  above 

3  P 
has  resulted,  since  the  assumption  that  —  =  0 

3Z 
prevents  the  operation  of  an  axial  pressure  gradient 

which  can  act  as  a  driving  force  for  the  z  component 

of  the  velocity. 
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ANALOG  SYMBOLS 


-(k3^Xj^(t)  +  k2X2(t) 


X2(t) 


Summing  Amplifier  and  Constant  Multiplier 


X(t) 


ki  X(t) 


Scale-factor  Potentiometer 


x(t)- 


Y(t^ 


X(t)  Y(t) 


Multiplying  Device 
106 


107 


U 

H 

+ 


^^o 

+J 

^-^ 

CN 

X 

.a^ 

CM 

n 

M 

c 

0 

+ 

•H 

■p 

^"^ 

•^ 

+J 

•0 

^-' 

G 

H 

0 

X 

u 

H 

;!>< 

c 

H 


o 


O 
•H 


I 
Id 

M 

;3 


X 


CN 

X 


APPENDIX   D 


PHYSICAL  PROPERTIES 

This  appendix  contains  the  numerical  values 
of  the  physical  properties  necessary  for  the  evalua- 
tion of  the  equations  presented  in  the  body  of  this 
dissertation.   Table  4  (JL9,  p.  535)  and  Table  5 
contain  some  of  the  properties  necessary  for  the 
solution  of  the  problem.   From  the  listings  in  Table 
4,  Figures  26,  27,  and  28  were  constructed  showing 
the  relationship  of  specific  volume,  viscosity,  and 
thermal  conductivity  of  steam  as  a  function  of 
temperature.   In  addition,  these  figures  display  some 
of  the  numerical  values  of  the  constamts  necessary 
for  substitution  into  Equations  (82),  (83),  and  (136). 
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TABLE  4 
PHYSICAL  PROPERTIES  OF  STEAM  AT  ATMOSPHERIC  PRESSURE 


T 

P 

U 

k 

Pr 

(F) 

(lb^/ft3) 

(ib^/ft 

sec> 

(BTU/hr 
ft  F) 

212 

0.0372 

0.870 

0.0145 

0.96 

300 

0.0328 

1.000 

0.0171 

0.95 

400 

0.0288 

1.130 

0.0200 

0.94 

500 

0.0258 

1.265 

0.0228 

0.94 

600 

0.0233 

1.420 

0.0257 

0.94 

700 

0.0213 

1.555 

0.0288 

0.93 

800 

0.0196 

1.700 

0.0321 

0.92 

900 

0.0181 

1.810 

0.0355 

0.91 

1000 

0.0169 

1.920 

0.0388 

0.91 

1200 

0.0149 

2.140 

0.0457 

0.88 

1400 

0.0133 

2.360 

0.0530 

0.87 

1600 

0.0120 

2.580 

0.0610 

0.87 

1800 

0.0109 

2.810 

0.0680 

0.87 

2000 

0.0100 

3.030 

0.0760 

0.86 

o 


g. 


Ill 


TABLE    5 
PHYSICAL   PROPERTIES 


Quantity  Value  Units  Refer- 

ence 


h^  970.3  BTU/lbjjj  (23) 


59.83  Ibjj/cu   ft  (19) 


T  212.  F  (28) 

sat 


32.1739  (Ibjn  ft/lbf   sec^)  (19) 


N  58.9    (lOOC)  dynes/cm  (29) 


B  1.2x10"^  1/C  (29) 

V 
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EXPERIMENTAL  DATA 

TABLE  6 
DROPLET  VAPORIZATION  TIMES 


Initial 

Vaporization 

Volume 

Time 

(cc) 

(sec) 

Initial      Vaporization 
Volume  Time 

(cc)  (sec) 


Average  Plate  Temperature  -  602  F 
Plate  Surface  -  304ss.  flat,  satin  finish 

1.3  417  0.52  294 
2^9  547  1.0  374 

4.4  636  2.6  535 
1.3  430  3.3  586 
3.7  615  4.2  602 
0.36  264  0.045  101 
1.9  480  1.6  479 
0.72  331  2.25  522 
0.18  181  2.0  493 
0.09  140  2.9  567 
0.19  193  3.75  578 
0.31  236  0.055  115 
0.60  321  0.04  89 
0.97  400  0.92  390 


Average  Plate  Temperature  -  605  P 

Plate  Surface  -  Graphite,  flat,  satin  finish 

0.44           250           1.5  403 

0.95           370           3.5  577 

0.51           280           2.4  499 

0.80           342           1.75  453 

1.1            375           4.5  602 

2.3  550  
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TABLE  6   (Cont'd) 


Initial 
Volume 
(cc) 


Vaporization 
Time 
(sec) 


Initial 
Volume 
(cc) 


Vaporization 
Time 


Average  Plate  Temperature  -  608  F 
Plate  Surface  -  304ss,  1  deg.  apex  angle,  satin  finish 


395 
380 
337 
303 
140 
100 
102 
282 


0.33 

229 

0.98 

0.45 

286 

0.87 

0.90 

386 

0.72 

0.10 

147 

0.61 

0.60 

314 

0.11 

0.20 

205 

0.045 

0.05 

100 

0.05 

0.83 

345 

0.42 

Average  Plate  Temperature  -  1014  F 
Plate  Surface  -  304ss,  1  deg.  apex  angle,  satin  finish 


172 

192 
203 
133 
123 
114 
155 
154 
178 
200 


0.39 

143 

0.59 

0.96 

213 

0.80 

0.81 

198 

0.97 

0.32 

139 

0.33 

0.60 

186 

0.26 

0.14 

96 

0.20 

0.08 

75 

0.50 

0.06 

61 

0.49 

0.05 

57 

0.69 

0.04 

55 

0.89 

APPENDIX   P 


DROPLET  SHAPE  UNDER  VARIABLE  GRAVITATION 
The  shape  of  a  liquid  droplet  with  a  surface 
tension,  N,  in  a  gravitational  field.  ^  ,    can  be 
found  by  solution  of  the  equilibrium  equation  (27,  p. 
444) 

1       1      AP 
ri     r2      N 

for  the  droplet  shown  in  Figure  29.   The  notation 
used  in  reference  (27)  is  altered  slightly  for  conven- 
ience in  the  problem  at  hand.   By  expressing  r^  and 
r^    in  terms  of  p  and  introducing  the  new  variables 

J   =  sin  e  ,  ^g^j 

*^  =  ^D  ^  =^   '  (162) 

the  equilibrium  equation  takes  on  the  forms 


dx     X       N 


dz 
dx 


(1  -  j2)^ 
119 


(163) 


(164) 
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Fig.  29.   Schematic  of  Water  Droplet. 
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Equations  (163)  and  (164)  are  integrated 
numerically  starting  from  the  upper  point  A  shown  in 
Figure  29.   At  this  point,  r^^  =  r^   from  symmetry, 
and  Equation  (160)  takes  on  the  form 

r^  -  .  (165) 


Pd 


r  2 


For  a  given  value  of  the  surface  tension,  N,  the 
shape  of  the  droplet  is  now  solely  determined  by  the 
pressure  at  point  A  which  is  directly  related  to  the 
radius  at  point  A  by  the  above  equation.   The  pres- 
sure and  radius  at  point  A  are  set  by  arbitrarily 
letting  z  =  d  at  point  A.   Therefore, 


r^  =    ^^    .  (166) 


It   d  is  chosen  small,  r,  becomes  large  and 
the  droplet  taXes  on  a  pancake  or  flat  spheroidal 
shape.   If  d  is  chosen  large,  r^  becomes  small  and 
the  droplet  takes  on  a  spherical  shape.   By  choosing 
numerical  values  of  d  cmd  solving  each  problem,  the 
shape  of  the  droplet  for  various  volumes  is  determined, 

Further  simplifications  to  the  equations 
occurs  if 
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E 


2   _ 


N 


(167) 


Starting  with  an  assumed  value  of  d,  the  integration 
is  begun  by  making  the  first  element  of  the  meridian 

At  the  end  of  this  arc,  see 
Figure  30, 


curve   r ■,  A  0  =  x 


dfl  + 


d  A  X 


(Ax) 


4E' 


(168) 


(169) 


2E'' 


After  the  values  of  z,  and  j,  are  found  from 
the  above  equations,  as  shown  in  Figure  30,  the  com- 
plete solution  is  marched  out  by  means  of  Equations 
(163)  and  (164)  as  follows: 


i+ 


1  -  Ji 


dz 

dj 

dx 

dx 

i 

dz 
dx 

AX 


(170) 


'i+1 


=  z.   + 


dz 
dx 


AX 


(171) 
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Fig.  30.   Path  of  Numerical  Integration. 
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where  Equations  (163)  and  (164)  take  on  the  forms 


dj_ 
dx 


■>i. 
^i 


X172) 


dz 
dx 


'       a  -  i,^)"^ 


(o^e  "^90) 


(173) 


dz 
dx 


-h 


(1  -  Ji^)^ 


(90«=9^130>.    (174) 


Note  that  in  this  numerical  integration  ©  =90  results 
in  a  calculational  singularity;  consequently,  this 
point  is  excluded  from  the  integration.   The  position 
X  is  defined  by 

dz 
dx 


'i+1 


=   X.   + 


AX 


(175) 


dz 
dx 


The  calculations  are  continued  up  to  the  point  B, 
where  the  meridian  cuirve  has  the  horizontal  tangent, 
BC. 
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The  volume  of  the  droplet  is  readily  computed 
during  the  integration  by  the  following  summation: 

n 

V  =   Z    ^^i  ■  2.  J  TT  (x.)2   .       (176) 
1=1  ^~-^  ^ 

The  average  droplet  thic)cness  is  calculated  from  the 
relationship 

V 
£  = •  (177) 

""    (x.     )2 
i,max 

The  preceding  equations  were  programed  on  the 
IBM  709.   The  results  are  shown  in  Figures  31  and  32. 
In  Figure  31,  the  calculated  values  of H  are  shown 
for  values  of  r  -  1  and  r  =  0.16.   Included  in  this 
figure  are  data  points  based  on  measured  values  of  V 
^"^  ^i.max  ^^0*°  reference  (5).   As  can  be  seen  from  this 
figure,  there  is  good  agreement  between  the  calculated 
and  measured  values  of  £  for  the  case  of  r  =  l.   Con- 
sequently, it  is  expected  that  the  calculated  values 
of  J.  for  ^   -  0.16  will  be  equally  accurate.   The  r  =  1 
curve  is  replotted  in  Figure  32  for  easier  determina- 
tion of  *■  at  lower  volumes. 
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Fig.  31.   Thickness  of  Water  Spheroid  as  a  Function  of 
its  Volume  for  r=l.o  euid  r=0.16. 
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Fig.  32.   Thickness  of  Water  Spheroid  as  a  Function 
of  its  Volume  for  r=l. 
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